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Abstract 

m ■ 

fvj ' Post-Newtonian expansions of the Brill-Lindquist and Misner-Lindquist solutions of the time- 

symmetric two-black-hole initial value problem are derived. The static Hamiltonians related to the 
expanded solutions, after identifying the bare masses in both solutions, are found to differ from 
^ ' each other at the third post-Newtonian approximation. By shifting the position variables of the 

Cn , black holes the post-Newtonian expansions of the three metrics can be made to coincide up to the 

0^ ■ fifth post-Newtonian order resulting in identical static Hamiltonians up the third post-Newtonian 

approximation. The calculations shed light on previously performed binary point-mass calculations 
at the third post-Newtonian approximation. 
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1 Introduction and Summary 



In a recent paper by the authors |]lj the claim was put forward that at the third post-Newtonian (3PN) 

approximation of general relativity point-mass models for binary systems have to be replaced by black- 

k> , hole models to become unique. To confirm the claim that binary point-mass models are incomplete at 

}^ • the 3PN approximation, in the present paper the Brill-Lindquist (BL) Q and Misner-Lindquist (ML) 

5^ I M,m solutions of the time-symmetric initial value problem for binary black holes are expanded into 

post-Newtonian series and post-Newtonian Hamiltonians related to these expansions are calculated. 

The BL and ML solutions are known to differ from each other topologically as well as geometrically 
P, 0. The interesting question therefore arises if at the 3PN order of approximation in the relative 
motion differences show up. The remarkable outcome of our calculations from Sec. 2 is that the two 
solutions have different Hamiltonians starting at the 3PN order, but also that these Hamiltonians can be 
made to coincide by shifting the centers of the black holes. Especially interesting shifts are those where 
the two black holes of the ML solution obtain vanishing dipole moments as to possibly coincide with the 
monopolar black hole potentials of the BL solution to higher orders. 

It is evident that at the 3PN order of approximation point-mass models in many-body systems are 
no longer applicable. In Sec. 3 we perform static (i.e. with linear momenta of the bodies and transverse- 
traceless part of the three metric set equal to zero) binary point-mass calculations using different reg- 
ularization methods which lead to different metric coefficients and Hamiltonians, i.e. we end up with a 
static point-mass ambiguity. Sec. 4 is devoted to some consistency calculations for the regularization 
procedures of Sec. 3 and to the comparison with previous results. Some details of the calculations are 
given in two appendices. 
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Although it holds that through the postulate of vanishing black-hole dipole moments in isotropic 
coordinates a unique static Hamiltonian can be obtained at the 3PN order of approximation, the problem 
of finding a similar unique total Hamiltonian is more complicated because the point-mass ambiguity 
detected in the previous paper ||l| is a dynamical one which includes the momenta of the objects as 
well as the radiation degrees of freedom of the gravitational field (transverse-traceless part of the three 
metric) . This ambiguity lies far beyond any known exact or approximate solutions of the Einstein field 
equations and thus, it can not be resolved in a way the BL and ML solutions for the constraint equations 
allow for a clear identification of the ambiguous static contributions. The static ambiguity has not been 
mentioned in the paper ||I| . 

We use units in which 16ttG = c = 1, where G is the Newtonian gravitational constant and c the 
velocity of light. We employ the following notation: x — (a:') {i = 1,2,3) denotes a point in the 3- 
dimensional Euclidean space endowed with a standard Euclidean metric and a scalar product (denoted 
by a dot). Letters a and b are body labels {a,b = 1,2), so x^ denotes the position of the ath body, and 
rua denotes its mass parameter. We also define Va :— x — Xa, ra := |ra|, no := Va/ra', and for a ^ b, 
Tab '■= Xa — Xb, Tab '■= \TCab\, ^^ab '■— ^ab/^ab] \ ' | stauds here for the Euchdean length of a vector. Indices 
with round brackets, like in 4'{2)i give the order of the object in inverse powers of the velocity of light, in 
this case, 1/c^. We abbreviate (5 (x — x^) by 5a- 

2 PN expanded Brill-Lindquist and Misner-Lindquist solutions 

In the static case, as defined above, the three metric can be put into conformally flat form [cf. Eq. (3) in 

0] 

m= (l + ^-^j ^^r (1) 

We use the representation of the BL and ML solutions as given in Ref. Q. For the BL solution the 
function cj) from Eq. (|l|) equals [cf. Eq. (3) in 1^] 

0^^=8f^ + ^V (2) 

Here the positive parameters ai and a2 can be expressed in terms of the bare masses mi and 7712 of the 
black holes and the coordinate distance ri2 between them. The relations are given in Appendix A. The 
ML solution is described by the function [cf. Eq. (4) in jgl] 

where Xi in ri = |x — xi| is the position of the center of the black hole 1 of radius a and X2 in r2 — |x — X2I 
is the position of the center of the black hole 2 of radius 6, relative to a given origin in the flat space; d„ 
(n > 2) are the positions of the image poles of black hole 1, e„ (n > 2) are the positions of the image 
poles of black hole 2, a„ and 6„ {n > 2) are the corresponding weights. For the ML solution the choice 
of the bare masses is not as obvious as in the case of the BL solution Q . We use the definition of the 
bare masses introduced by Lindquist in Q. Then the solution (^ can be iteratively expressed in terms 
of the bare masses mi, 1712 and the vector ri2 connecting the centers of the black holes, cf. Appendix B. 
In the following we identify the bare masses of both solutions. 

The post- Newtonian expansions of the functions cj)^^ and (j)^^ can be written as follows 

^^^ = ^'f-ri + (t^fi) + '/'fe) + ^f&) + (t>fio) + ^ ( ^ j ' (4) 



where 0(„) are functions of ^PN order, as they belong to the three metric. The details of the expansions 
are given in Appendices A and B. The results read 
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The equations (|q)-(|l^) show that the ML solution at the 4PN order of approximation attributes 
a dipole moment to each black hole whereas the BL solution, as it is already evident from the exact 
expression (0) , shows monopoles only. In shifting the centers of the black holes one can arrange that also 
in case of the ML solution the dipole moments do vanish. To show this let us introduce 
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Then the shifted ML solution can be defined as 

^(x;xi,X2) := (/)<|i.N(x;xi +ari2,X2 +/5r2i) 
where a and (3 are some dimensionless parameters. We have found that for 
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the shifted ML solution coincides with the BL solution up to the 5PN order of approximation: 



I ML shifted/ , \ /BL / , n 

»<5PN (x;Xi,X2) = 0<5pn(x;Xi,X2) 



The Hamiltonian we calculate by means of formula 

H = — lim i) dai 
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where S (0, R) is a sphere of radius R centered at the origin of the coordinate system. Making use of Eqs. 
(g) and (^) we obtain that the BL and ML solutions lead to the Hamiltonians 

ifBL ^ 327r(ai+a2), (20) 

OO 

i/ML ^ 327r(a + 6) + 327r^(a„ + fe„). (21) 

n=2 

Using Eqs. ( pO| ) and (|T]) [or Eqs. (|q)-(p^ together with Eq. ([l9|)] we calculate the static Hamiltonian 
up to the 3PN order of approximation (notice: the nPN Hamiltonian is determined by the {n + 2)PN 
three metric). We obtain, dropping the total mass rrii + m2 contribution (in the reduced variables 0) 

i^lsPN - - jl^Jda, (0|U), = -; + 2^^ - i(l + '^);i + ^(1 + 3^)^' (22) 

S(0,i?) 

ff^g'^P^ = - ^lim ^da, (</,^U), = -; + 2^ - i(l + -);i + ^(1 + 2-)^- (23) 

Obviously, the both 3PN Hamiltonians are different. The difference vanishes, however, if in the case of the 
ML solution the shifted solution for the potential function is used in the calculation of the Hamiltonian. 
Therefore, the postulate of vanishing black-hole dipole moments in isotropic coordinates yields a unique 
static 3PN Hamiltonian for binary black holes. The condition of vanishing dipole moments needs the 
metric coefficients; it can not be formulated on the Hamiltonian level alone. 

3 Binary point-mass calculations 

In this section we show the results of the static binary point-mass calculations. In the static case the 
Hamiltonian constraint equation for the two-body point-mass system in the canonical formalism of ADM 
reads 

(l + ^A^cl>=-Y,T^aSa (24) 

^ ^ a 

(for binary systems all sums run over a = 1,2). Eq. (p3) yields the following formal expansion 

^ ^ a \ a / ri=0 ^ ^ 

Using Eq. ( P5|) we obtain the Hamiltonian constraint equations valid at individual orders in 1/c. They 
read 

A0(2) = ~^maSa, (26) 

a 

A0(4) = ■^(f>{2)^maSa, (27) 

a 

A<^(6) = (-^'^(2) + ^<^(4))E^"-'^- (28) 

^ ^ a 

^^'^^'> " (^'?^?2) - ^'/'(2)'/'(4) + g0(6) j Xl"^"'^'" (29) 

^ ^ a 



All Poisson equations (^6|)-(30) are of the form 



(31) 



where the function / is usually singular at x = x^. One can propose three different ways of solving 
equations of type (p|). The first two ways are based on the following sequence of equalities: 



A-l (^/(x)5„Ua-1 (^/,.eg(x,)U =^/,,g(x„)A-Ma = -i^^/,.eg(x,)l 
\ a / \ a /a a 



(32) 



where /rog(xa) is the regularized value of the function / at x = x^, defined by means of the Hadamard's 
"partie finie" procedure ||]. The difference between the two first methods relies on different evaluating 
the regular value of the products of singular functions. In the first method we use 



(/l(x)/2(x))(5a = /lrog(Xa)/2rcg(Xa)(5a, 

whereas in the second method instead of the rule (gSJ) we apply 



(33) 



(34) 



For the developments in the book by Infeld and Plebahski [^ , it was crucial that the both regularization 
procedures coincided, i.e., (/i/2)rcg = /ireg/2rcg ("tweedling of products"). In the third method we 
regularize the Poisson integral rather than the source function: 
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Let us denote the results of applying the regularization method based on Eqs. (p2) and (p3) by primes, 
the results of the method based on Eqs. ( p^ ) and ( p4| ) by double primes, and the results of the method 
based on Eq. ( p^ ) by triple primes. It turns out that up to 1/c^ the results of the three methods coincide 
and, moreover, they are identical with the results of the expansions of the BL and ML solutions: 
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The results of the first method coincide with the expansion of the BL solution, what we have checked up 
to the 5PN order: 

^' --^fL, 71 = 2,4,6,8,10. (37) 
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reads 



*('lO) 



,BL 
»(10) 



1 



mTmX 



2(167r)5 



mi 
ri 



7712 

r2 



The functions 



'(8) 



and 



'(10) 



are equal to 



'(8) 



(and 



i?g'^), and the 



(38) 



"(8) 



iiii 

3(10) 



/BL 



,BL 

»(10) 



1 



raxTfii 



2(167r)4 
1 



''12 

mim2 



2 "1 • "12 2 "2 • ni2 

nil 2 "^2 2 



"2 -1112 



4(167r)5 



'12 



m2{m2 + 677T,i) 2 m^{mi + 677^2) 



ni • ni2 



rr 



(39) 



?71l77l2 

ri2 



(2mi + 1712) (- (2?7i2 + rni) — 

n r2 



(40) 



Notice, the mass parameters 7711 and 7712 in Eqs. ( p6)-(|4C| ) denote total rest masses of infinitely separated 
bodies whereas the mass parameters in Eqs. (P4|)-(30[) denote some formal rest masses only; the former 



results from the latter by our regularization procedures. The total rest-mass parameters we identify with 
the bare masses of the BL and ML solutions. 

By means of Eq. (|l9| ) we calculate the static Hamiltonian up to the 3PN order applying the functions 

»'<^('„)' and </)'(';). 



^' '^" and (j)'/^)- The results are (dropping the total mass mi + m2 contribution): 



H'SI^^' = - ^^Jd<y^ (^'<5Pn),, = -^ + ^ - i(l + -)^ + ^(1 + 3-)^: (41) 

S(0,R) 

H'^i^r = l^Jda.. (<^^',pn)^^ = -^ + i - i(l + ^);i + ^(1 + h^^- (43) 

S{0,R) 

Obviously, at the 3PN order of approximation, the three Hamiltonians differ from each other and from 
the Hamiltonian ([2^) obtained from the expanded ML solution. The Hamiltonian ( [4l| ) coincides with the 
Hamiltonian ( |2^ ) obtained using the PN expansion of the BL solution, what obviously follows from Eq. 

(§)• 

4 Consistency calculations and comparison with the previous 
results 

In the region Q := B (0, R) \ [B (xi, £1) U B (x2, £2)] (where B (xa, Sa) (a — 1, 2) is a ball of radius Sa 
around the position Xq of the ath body and B (0, R) is a ball of radius R centered at the origin of the 
coordinate system) the right-hand sides of the Eqs. (p6|)-(po|) vanish, so the functions </)(„) fulfil the 
Laplace equation in this region. Applying Gauss's theorem we thus obtain 

=: / d^o; A(/)(„) = j> da.i (j)(ri),i + j> da^ 0(„)^i + j> dat (/)(„),,;, (44) 

n as(xi,£i) aB(x2,£2) dB(0,R) 

with the normal vectors pointing inwards the spheres dB (xa,£a) and outwards the sphere dB(0,R). 
From Eq. (^) it follows that 

lim j) dai(j)(^r^).i = - Van i dcTi (/)(„) ^^ - lini ido-j0(„)_i, (45) 

dB{0,R) as(xi,ei) 9B(x2,e2) 

SO the Hamiltonian i?„pN at the ?iPN order, according to Eq. (|l9|), can be calulated as 

HnPN = - lim i dai (f>(2n+'k)A ~ 1™ f doi 4>(2n+i),i- (46) 

ei^O J £2^0 J 

9-B(xi,ei) 9B(x2,e2) 

We have used Eq. (46) and the functions (j)', x, 0'/ s, and 0'/'-,, to calculate the static Hamiltonian up to 



^(n)' ^(n)' "^"^ ''^(n)' 



the 3PN order. The integrals over the spheres dB (:x.a,£a) diverge as Sa — > 0, so to calculate them we 



have used the Hadamard's procedure M. The results coincide with those given by Eqs. (41)-(|4 
The nPN Hamiltonian can also be written in the form of a volume integral 



HnPN = d^x/^<j)(2n+A), (47) 

so still another way of calculating it relies on direct integration of the (minus) right-hand sides of Eqs. 
(p6|)-(p0|). To do this one must use the Hadamard's regularization [0 together with the rule (|3^) or (p3). 



For the functions 0'/ s we have used the rule ( p3| ) (as in deriving the functions </'f„')), whereas for the 



functions (/ )"„^ and (j)'!^^ we have apphed the rule (|4|). The results coincide again with those given by 
Eqs. (ill)-®. 

In the paper |l| we applied the 'double prime' regularization procedure and had thus obtained the 
Hamiltonian (|4|); see Sec. VI in Ref. [Q. The same result we had obtained from the static n-body 
Hamiltonian of Ref. |L0| by applying some expansion-and- limiting procedure. If we take the n-body 
static Hamiltonian of |10[ but specialize it to the two-body case adapting only those terms which are 
directly finite, we get the BL result of Eq. ([2^). It corresponds to our 'single prime' regularization 
procedure described above. 
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A Brill-Lindquist solution 

We use the representation of the Brill-Lindquist (BL) solution taken from Appendix A of g. The BL 
solution can be written in the form [cf. Eq. (3) in [|6|] 

^BL ^ 8 (^^ + £^y (48) 

The bare masses mi and 7712 of the black holes depend on the parameters ai, a2 and the coordinate 
distance ri2 between the black holes [cf. Eqs. (A7) and (A8) in [^]: 

mi = 327rai M -h — ) , TO2 = 327ra2 ( 1 + — ) . (49) 

The unique positive solutions of Eqs. (^9|) for ai and q;2 read 

«2 - -J(^2ri2 + ^^ij^)+^ri2l/4+^:i— ^+(::^ij-^ (51) 

We expand the right-hand sides of Eqs. (^y) and (|5l| ) in powers of 1/c taking into account that the 
masses mi and m2 can be regarded as being of order 1/c^. The results we substitute into Eq. ( [l8|) to obtain 
the post-Newtonian expansion of the function (f)^^ . Such obtained functions 0?^^ for n = 2,4, 6, 8, 10 are 
given in Eqs. @, §), @, (|), and (0), respectively. 

B Misner-Lindquist solution 

The form of the Misner-Lindquist (ML) solution we use is taken from Appendix B of Q . The ML solution 
we write in the form [cf. Eq. (4) in jg]] 

^--(^^)-|(^^^). 

where Tq — |x — Xa | {a — 1,2) and Xa is the position of the center of the ath black hole relative to a given 
origin in the flat space, a and b are the radii of the black hole 1 and 2, respectively; d„ {n > 2) are the 
positions of the image poles of black hole 1, e„ {n > 2) are the positions of the image poles of black hole 
2, a„ and &„ (n > 2) are the corresponding weights. 
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The black hole 1 together with its odd images and the even images of the black hole 2 are located 
on the positive z axis, whereas the black hole 2 together with its odd images and the even images of the 
black hole 1 lie on the negative z axis. The relative distances |x — d„| and |x — e„| entering Eq. ( p2| ) 
can be expressed by radii a and h of the black holes and their relative position vector X\2 as follows (here 
n > 2) 

|x - d„| 
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for n odd, 
for n even, 
for n odd, 
for n even. 
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The quantity /xq entering Eqs. (54) is given by 
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The weights a„ and 6„ (n > 2) of the image poles depend on the radii a and b of the black holes and the 
distance ri2 between them: 
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for n even, 
for n odd, 
for n even, 
for n odd. 



(56) 



To obtain the post-Newtonian expansion of the ML solution we expand the right-hand side of Eq. 
(52) in powers of 1/c taking into account that the radii a and b are of order l/c^. To do this we use Eqs. 
(53)-(p6[). We have found that only the first five terms from the sum on the right-hand side of Eq. ( p2|) 
contributes to the 5PN order. The result of the expansion reads 
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For the ML solution we use the definition of the bare masses mi and ?7i2 introduced by Lindquist in 
, The masses rrii and TO2 depend on the radii a and b of the black holes and the distance ri2 between 



the centers of the black holes [cf. Eqs. (B19) and (B20)]: 
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smh 2/io 7 n < . j — — 7 H — > , (59) 

ri2 ^ |_sinh2T7,^o smh[2(n + l)//o — ^^2] sinh[2(?i — Ij^o — M2J J 

where fj,2 is given by 

sinh/i2 = — sinh2/io- (60) 

?'12 

We have iteratively solved Eqs. (Bq) and (bm with respect to a and &. The result, valid up to the 5PN 
order of approximation, read 

1 nil 1 TO1TO2 1 mi?7i2(2mi + 37712) 1 mim2{mi + Amim2 + irn^) 



16tt 2 (167r)2 2ri2 (167r)3 8r^2 (IStt)" 8r? 
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1 mim2{2m,l + 14rnlrn2 + ISmimj + 5rnl) ^n f ^ \ (c^^ 



(167r)5 32rt2 

the equation for b can be obtained from the above one by replacements mi —^ ?7i2 and 1712 —^ mi. 

To obtain the PN expansion of the ML solution we substitute Eqs. ( |6l|) into Eq. (|57|). Such obtained 
functions </«^^ foj. „ ^ 2, 4, 6, 8, 10 are given in Eqs. (§), (0), (|), (|ll|)' and (|l|), respectively. 
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